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We calculate the axial couphngs of mesons and baryons containing a heavy quark in the static 
limit using lattice QCD. These couplings determine the leading interactions in heavy hadron chiral 
perturbation theory and are central quantities in heavy quark physics, as they control strong decay 
widths and the light-quark mass dependence of heavy hadron observables. Our analysis makes 
use of lattice data at six different pion masses, 227 MeV < < 352 MeV, two lattice spacings, 
a = 0.085, 0.112 fm, and a volume of (2.7 fm)^. Our results for the axial couplings are gi — 0.449(51), 
g2 ~ 0.84(20), and gs = 0.71(13), where gi governs the interaction between heavy-light mesons and 
pions and 512, 3 are similar couplings between heavy-light baryons and pions. Using our lattice result 
for g3, and constraining l/mg corrections in the strong decay widths with experimental data for Ec*^ 
decays, we obtain Tp^'^^Ab 7r±] = 4.2(1.0), 4.8(1.1), 7.3(1.6), 7.8(1.8) MeV for the E+, E", E*+, 
initial states, respectively. We also derive upper bounds on the widths of the H^'*'' baryons. 



Introduction. — Significant progress has been made in 
the last few years in uncovering the spectrum and de- 
cays of hadrons containing heavy quarks at the dedi- 
cated B factories, the Tevatron, and the LHC. Accurate 
lattice QCD calculations are required to confront data 
from these experiments with the Standard Model. These 
lattice calculations involve extrapolations in the masses 
of the light quarks, which require theoretical guidance. 
For hadrons containing a single heavy quark, the rel- 
evant effective theory is known as heavy-hadron chiral 
perturbation theory (HHxPT) [iHSI , which is built upon 
two of the most important symmetries of QCD: chiral 
symmetry and heavy-quark symmetry. At leading or- 
der, the HHxPT Lagrangian contains three axial cou- 
plings gi, g2, and 53. The coupling gi determines the 
strength of the interaction between heavy-light mesons 
and pions, while 172 and 33 similarly determine the inter- 
action of heavy-light baryons with pions. Consequently, 
these couplings are central to the low-energy dynamics 
of heavy-light hadrons, and can be used to calculate the 
widths of strong decays such as SJ^ ' — ^ A^ tt. The ax- 
ial couplings are calculable from the underlying theory 
of QCD, using a lattice regularization. The mesonic cou- 
pling gi has been previously studied in lattice QCD with 
Nf = or Nf = 2 dynamical quark flavors [MS]. In the 
following, we present the first complete calculation of gi , 
g2, and 33 in Nf = 2 + 1 lattice QCD, controlling all 
systematic uncertainties. We use our results to calculate 
r[I][*-' Ab TT^] and give bounds on r[S'^'*'' Sf, tt]. 
Technical details of the analysis that are omitted here 
for brevity will be presented in a forthcoming paper. 

Lattice QCD calculation. — The heavy hadrons consid- 
ered in the lattice calculation are the lowest-lying heavy- 
light mesons and baryons containing light valence quarks 
of the flavors u or d. We work in the heavy-quark limit 
niQ = 00 where the axial couplings are defined, and as- 



sume isospin symmetry. The heavy-light mesons occur in 
degenerate pseudoscalar and vector multiplets, described 
by interpolating fields ~ and P*^ ~ 07ai9% 

where is a light quark of flavor i and Q is a static 
heavy antiquark. In the heavy-light baryon sector, we 
include both the states with si — and s; = 1, where 
si is the (conserved) spin of the light degrees of free- 
dom. The states with si = I are described by an in- 
terpolating field S^^^ Eabc {Cjf,)fsj qlfs ql^ Qca that 
couples to the isotriplet states with both J — 1/2 
(Eg) and J = 3/2 (Sg), which are degenerate in the 
heavy-quark limit. The isosinglet si — baryon Aq 
has J — 1/2 and is described by an interpolating field 
^ f-abc (C'75)/37 9a/3 ^l-y Qca- The axial couplings can 
be extracted by calculating matrix elements of the axial 
current ~ dj^jbu: 

{P^\A,\P,)=-2{g,Ue;, 

{Sdd\A^,\Td^:) ^-{g:iUU^,U. (1) 

Here, v is the four-velocity, is the polarization vec- 
tor of the P* meson, U is the Dirac spinor of the T 
baryon, and the U*^ 's are the "superfield spinors" of the 
S baryons TO . At leading order in the chiral expan- 
sion, the "effective axial couplings" {gi)cS defined via 
([1]) are equal to the axial couplings gi that appear in 
the HH^PT Lagrangian. The next-to-leading-order ex- 
pressions for {gi)ca are given in Ref. pUj. To calculate 
the matrix elements ([T]) in lattice QCD, we set v = 
and construct Euclidean two- and three-point correla- 
tors CH{t) = (xH(x,i) xL(x,0)) and Cn^H'it,!') = 
E^'iXH'M A^{^',f) xl,(x,0)), where t > i' > and 
Xh are the interpolating fields of the heavy hadrons as 
defined above. We form the ratios 

RiM^- ' (2) 
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R2{t,t') = 2 



(3) 



and the double ratio (needed because of the nonzero S—T 
mass sphttmg) 



(4) 



Here, ^, p are the Lorentz indices from the axial current 
or the interpolating fields for P* and S and are summed 
over when repeated. Using ([T]) and the spectral decom- 
position of the correlators, one finds that 



0{e 



(5) 



where the 5i are related to the energy gaps of the lowest 
contributing excited states. 

The calculations presented in this work make use 
of lattice gauge field configurations generated by the 
RBC/UKQCD collaboration [H] with 2 + 1 flavors of 
light quarks, implemented with a domain- wall action that 
realizes lattice chiral symmetry. The details of the en- 
sembles included in our analysis can be found in Ta- 
ble [ij We computed domain- wall light-quark propagators 
for a range of unitary (am. 



(val) 



d'rn'u'^f') and partially 



quenched {arn^^'^^^ < ^'m'^'^f') quark masses. As shown in 
the lower part of the table, we have data with (valence) 
pion masses ranging from 227 to 352 MeV, two lattice 
spacings, a — 0.085, 0.112 fm, and a large lattice volume 
of (2.7 fm)'^. The sea-strange-quark masses are about 
10% above the physical value, and we assign a 1.5% sys- 
tematic uncertainty to our final results to account for 
this, based on the size of the effect on similar observables 
as studied in Ref. [11]. For the light-quark propagators, 
we used gauge- invariant Gaussian smeared sources to im- 
prove the overlap of the hadron interpolating fields with 
the ground states. We constructed the three-point func- 
tions CH^H'(t,t') using light-quark propagators with 
smeared sources at (x, 0) and (x, t) and a local sink at the 
current insertion point (x', t'), for various separations t as 
shown in Table |T] The bare lattice axial current requires 
a finite renormalization Za to match the continuum cur- 
rent, = Za ^7^75^. We used nonperturbative results 
for Za obtained by the RBC/UKQCD collaboration [TT] . 

The action for the static heavy quark is a modified 
form of the Eichten-Hill action [12] in which the standard 
gauge links are replaced by HYP (hypercubic) smeared 
[13] gauge links, resulting in improved statistical sig- 
nals for the correlators [Tl]. To study heavy-quark dis- 
cretization effects and optimize the signals, we generated 
data for tihyp = 1,2,3,5,10 levels of HYP smearing, 
corresponding to different lattice actions for the heavy 
quarks. These actions have the same continuum limit, 
but may scale differently. Our final analysis focuses on 
^^HYP = 1, 2, 3. 



Ensemble 


a (fm) 


X T 


(sea) 


ml"=' (MeV) 


A 


0.1119(17) 


24^ X 64 


0.005 


336(5) 


B 


0.0849(12) 


32 X 64 


0.004 


295(4) 


C 


0.0848(17) 


32 X 64 


0.006 


352(7) 


Ensemble 


(val) 


^r'> (MeV) 


mi""' (MeV) 


t/a 


A 


0.001 


294(5) 


245(4) 


4, 5, 10 


A 


0.002 


304(5) 


270(4) 


4, 5, 10 


A 


0.005 


336(5) 


336(5) 


4, 5, 10 


B 


0.002 


263(4) 


227(3) 


6, 9, 12 


B 


0.004 


295(4) 


295(4) 


6, 9, 12 


C 


0.006 


352(7) 


352(7) 


13 



TABLE I. Details of gauge field ensembles (upper section, 
see also Ref. [Til!) and "measurements" (lower section). The 
superscripts v, s on indicate the masses of the quarks in 



the pions, equal to am^J^^ or am'^^ 
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FIG. 1. Ratios Ri{t,t') as a function of the current insertion 
time slice t' , for t/a = 10, at a = 0.112 fm, am^"^'' = 0.002, 
n.HYP — 3. 



In Fig.[T] we show examples of numerical results for the 
ratios ([2]), ([3|, and We observed plateaus in Ri{t, t') 
as a function of t' , and we averaged the ratios in this re- 
gion, which is essentially equivalent to taking Ri{t,t/2). 
We denote these averages as Ri (t) . To obtain the ground- 
state contributions according to one needs to calcu- 
late limf_>.oo Ri{t). To this end, we performed fits of the 
data using the functional form Ri{t) = {gi)cff — Ai e~^^ * 
with parameters {g^)cf^^ Ai and (5i, depending on the lat- 



tice spacing a, the quark masses a'rn^J^\ '^''^u'^d'^ ^^'^ ^^e 
smearing parameter njjYp . This functional form only in- 
cludes the leading contributions from excited states, but 
was able to fit the data well, as shown in Fig. [2] We used 
the results and uncertainties for the gap parameters Si 
from the fits at the coarse lattice spacing to constrain the 
fits at the fine lattice spacing, where we have fewer val- 
ues of t/a. As explained in Ref. [TS], we then additionally 
constrained the parameters Ai (independently for the two 



FIG. 2. Fits of the t-dependence of R^{t), for a = 0.112 fm, 



(val) 



= 0.002, nHYP = 3. 
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different lattice spacings), using information from initial 
fits of data from ensembles A and B. This allowed us to 
perform fits using the same form of the function Ri{t) 
even for the data from ensemble C, where we have only 
one value of t/a. To estimate the systematic uncertain- 
ties caused by higher excited states, we calculated the 
shifts in {gi)cS at the coarse lattice spacing when remov- 
ing one or two data points with the smallest t/a (= 4, 5) 
or adding a second exponential to the fits ^15j . Repeated 
fits of Ri (t) for a bootstrap ensemble allowed the calcula- 
tion of the covariance matrices describing the correlations 
of the results for ((7i)eff from common ensembles of gauge 
field configurations. 

Having obtained the results for {gi)eS, we then per- 



(vv) 



formed fully correlated fits of the a-, mi- 
dependence. For {gi)cS, we used the function 



and m^-i^^^- 



(51)0 



.91 



I + /i(5i, mi--), m^-), L) + d,,nn.. 



„(vv) r^(vv)l2 , (vs) r (vs)l2 



(6) 



where 171, c"i^\ c"i^\ {di.nnYp} ^^'^ the free parameters. 
For (g2)cff and (53)0!?, we performed coupled fits using 



. (VV) r (vv).2 



'TT 

.(vs) 



m 



(vs)l2 



(7) 



(for i — 2,3), where the free fit parameters are 52, 93, 

4'"'^ 4'"'''' 4'"'^ 4''*'^ {c'2,«HYp:'^3,nHYp}- The functions 
fi in Q and ([7| are the nonanalytic loop contributions 
in partially quenched S'[/(4|2) HHxPT and can be found 
in Ref. [TD]. They also include the leading effects of the 
finite lattice size L (because of our large volume, the 
finite-volume corrections were smaller than 3% for all 
data points). The functions fi depend on the renormal- 
ization scale /i, but this dependence is canceled exactly 
by the /^-dependence of the counterterms c^^^'* and cf'^'^ . 



FIG. 3. The (real parts of the) fitted functions {gi)cs, (52)011, 
(.93)011, evaluated in infinite volume and uhyp = 3, for the 
unitary case m^^-' = ml^"' = m^. The dashed line corre- 
sponds to a = 0.112 fm, the dotted line to a = 0.085 fm, 
and the sohd line to the continuum Hmit. The shaded regions 
indicate the la statistical uncertainty. Also shown are the 
data points, shifted to infinite volume (circles: a — 0.112 fm, 
squares: a = 0.085 fm). The partially quenched data points 
(open symbols), which have m^^-' < m^''\ are included in 



the plot at m.^ 



,(vv) 



, even though the fit functions actu- 



ally have slightly different values for these points. 



The parameters rfj.nnYP fo^' ^ach tihyp describe the lead- 
ing effects of the nonzero lattice spacing, which are multi- 
plicative corrections proportional to as a consequence 
of the lattice chiral symmetry of the domain- wall action. 
In ([7]), the quantity A is the S —T mass splitting, which 
we set to A = 200 MeV in our fits, consistent with ex- 
periments [ini [H] and our lattice data (note that A does 
not vanish in the chiral or heavy-quark limits). 

To determine for which values of nuYP the order-a^ 
corrections in ^ and ^ adequately describe the lattice 
artefacts in the data, we started from fits that included 
all values of nnyp, and then successively removed the 
data with the largest values of tihyp- After excluding 
JT-HYP = 10 and riHYP = 5, we obtained good quality-of- 
fit values [Q = 0.70 for (gi)eff and Q = 0.92 for {92,3)cs], 
and the results were stable under further exclusions. Our 
final results for the axial couplings, taken from the fits 
with riHYP = 1, 2, 3, are 

51 = 0.449 ± 0.047 ,tat ±0.019 ^yst , 
52=0.84 ±0.20 Stat ±0.04syst, 
53 = 0.71 ±0.12 

stat 

±0.04 

syst • (8) 

Independent fits for each riHYP (1, 2, 3, 5, 10) gave re- 
sults consistent with p 
atic uncertainties in 



The estimates of the system- 
include the following [15^: ef- 



4 



fects of next-to-next-to-leading-order terms in the fits to 
the a- and mTr-dependence (3.6%, 2.8%, 3.7% for gi, g2, 
gs, respectively), effects from the unphysically large sea- 
strange-quark mass (1.5%), and effects from higher ex- 
cited states in the t ^ oo extrapolations of Ri{t) (1.7%, 
2.8%, 4.9%). The resulting mass- and lattice-spacing de- 
pendence of the effective couplings from the fits with ([6| 
and ([t]) is shown in Fig.[3j Note that the functions {g2)cS 
and ((73)cff develop small imaginary parts for pion masses 
below the S ^ Tn threshold at = A [TO] (the lat- 
tice data are all above this threshold) , and the real parts 
are shown in the figure. The fitted coefficients di, uhvp 
are consistent with zero within statistical uncertainties, 
and the analytic counterterms c|^^'' and cf^"^ are natural- 
sized (when evaluated at ^ = 47r/^ with — 132 MeV), 
indicating that the chiral expansions of {gi)cS are under 
control for the range of masses used here. 

Calculation of strong decay widths. — At leading order 
in the chiral expansion, the widths for the strong decays 
S' ^ T TT are 

nS^Tn]=c^-^(gs + ^y^\p.\^ (9) 

where S and T now denote physical s; = 1 and si = 
heavy baryon states such as Sf, and A;,, |p7r| is the mag- 
nitude of the pion momentum in the S rest frame, and 
Cf is a flavor factor, equal to 1 for Eg'' — > Aq tt*, 1/\/2 

for Sg*-* Sq tt*, and 1/2 for Sg*-* — > Sg tt'^. Here we 
modified the mg = oo expression for F |18) by includ- 
ing the term kj/uiq. Terms suppressed by (m^/A^)^ 
and (Aqcd/"^q)^, which are omitted from M, lead to 
small systematic uncertainties in F. To determine ki/2 
and K3/2, we performed fits of experimental data 19J for 
the widths of the Sj? (J = 1/2) and thc"E*++, 

E*° (J = 3/2) using M, where we constrained 33 to our 
lattice QCD result H and set mg = ^Mj/^. These 
fits gave Ki/2 = 0.55(21) GeV and K3/2 = 0.47(21) 
GeV. We then evaluated (joj) for toq = ^My to ob- 
tain predictions for the decays of bottom baryons. Our 
calculated widths F[E^*'' -> A;, tt*] as functions of the 

— Af, mass difference are shown as the curves in 
Fig. |4j Using the experimental values of the baryon 
masses HZl HU, our results for F[E^*^ ^ A^ TT*] in MeV 
are 4.2(1.0), 4.8(1.1), 7.3(1.6), 7.8(1.8) for the E+, E^, 
E^^, E^~ initial states, respectively, in agreement with 
the widths measured by the CDF collaboration ■ The 
decays E',^*^'^ E^n^, Sg tt" and E',^*^' ^ E^n+, 
may also be allowed, depending on the mass differences. 
With a spin-averag ed E[^*^ - Eb splitting of 153(21) MeV 
(based on lattice data from Ref. [5D]), and assuming 
M{El) - M{E[) « M(E*) - M(Eb) = 21(2) MeV [IS], 
we obtain upper bounds of 1.1 and 2.8 MeV (CL=90%) 
for the total widths of the SJ, and 5^, respectively. 

Conclusions. — We have presented a lattice QCD cal- 
culation of the axial couplings of hadrons containing a 
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FIG. 4. Widths of the decays E^*'=^ ^ TT* as functions of 
the E^*' - At mass difference. The curves (solid: E;,, dashed: 
Ej^) and shaded regions show our predictions and their uncer- 
tainties. The data points are from CDF [17| . 

heavy quark in the static limit, including for the first 
time the baryonic couplings. We have used these results 
to predict the strong decay widths of bottom baryons. 
Our calculation of the axial couplings controls all system- 
atic uncertainties by using two different lattice spacings, 
low pion masses, a large volume, and the correct next-to- 
leading-order expressions from HHxPT. Since the axial 
couplings are essential for chiral extrapolations of lattice 
data, their accurate determination is of broad significance 
in flavor physics phenomenology. 
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